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Abstract 

The light quark propagation in the confining vacuum, described by an (infinite) set 
of gauge-field vacuum correlators, is studied in detail. To keep gauge invariance at each 
step the system of light quark and a heavy antiquark is considered, and the nonlinear 
equations are written explicitly for the quark propagator in the limit of large Nc- For 
the lowest (Gaussian) correlator the system is studied in different approximations, and 
the relativistic WKB method is used to demonstrate the scalar confining interaction 
of light quarks, which implies chiral symmetry breaking. The chiral condensate is 
estimated by the WKB method, and connection to the density of global zero modes 
is clarified. The higher even order correlators are shown to yield the same properties 
of scalar confining interaction for light quarks. No attempt was made to solve the 
obtained nonlinear equations numerically, but the qualitative conclusion of connection 
between confinement and chiral symmetry breaking is drawn, and an estimate of the 
chiral condensate is performed. 



1 Introduction 



The light quark propagation in the QCD vacuum displays two phenomena: 
confinement and chiral symmetry breaking (CSB). The first one can be most 
clearly studied in the example of heavy quarks, where the linear confining 
potential (equivalent to the area law of the Wilson loop) is a good order 
parameter in the absence of dynamical quarks. Recently the new vacuum 
correlator method was introduced [1] which successfully describes dynamics 
of confinement, and the connection with the dual Meissner effect can be 
clearly seen (for a recent review see [2]). 

The confining quark dynamics for massive quarks was formulated in the 
vacuum correlator method, and the method allows to take into account spin 
degrees of freedom [3-5] but only as a perturbation in powers of 1/M. E.g. one 
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can extend the method to effectively calculate spin-dependent contributions 
to the mass of p meson, but one fails in the case of pion. 

This difficulty can be traced back to the lack of CSB in the Feynman- 
Schwinger representation for the quark propagator with spin [6]. 

To treat CSB one can use the most popular way - to consider a gas or 
liquid of topological charges in the vacuum. For the case of instantons it was 
argued that CSB occurs ( at least for large Nc) for any instanton density ([7], 
for a recent review see [8]), a similar result holds for dyonic gas [9]. 

Therefore a logical way of creating a natural QCD environments for light 
quarks in this approach consisted of placing topological charges inside a con- 
fining vacuum, described by vacuum correlators. It was shown in [10] that 
such a construction indeed provides CSB together with confinement in a 
gauge-invariant way, and moreover instantons are stabilized at large distances 
due to confinement [11]. In particular one could visualize in such vacuum the 
double nature of pion as a qq bound state and as a Nambu-Goldstone particle 
in one and the same physical Green's function. 

There are visible defects in such construction however. Firstly, it looks 
artificial to superimpose instantons in the confining background. Secondly, 
there is no explanation of why confinement and CSB occur together in the 
confining phase and why they disappear simultaneously above Tg, as it was 
shown repeatedly on the lattice [12]. 

There is another type of approach to CSB in connection with confinement 
[13], where Dyson-Schwinger equations are used for the isolated quark prop- 
agator with selfinteraction via modified gluon exchange However the system 
under study is not gauge invariant and physically there one does not take 
into account the QCD string, connecting the quark to antiquark. 

In this paper we choose therefore another way. We start with the QCD 
Lagrangian and derive from that the effective Lagrangian of light quarks, 
assuming that certain gluon field correlators are nonzero, which are known 
to yield confinement (i.e. linear potential) for static quarks. It is not clear 
from the beginning what will happen for light quarks ( with vanishing mass) 
and whether they would be confined at all. 

Our main concern in what follows is to keep Lorentz and gauge invari- 
ance. The effective quark interaction is nonlocal and to have gauge-invariant 
equation one should consider qq Green's functions. 

The simplest setting for which confinement and CSB can be studied in 
the gauge-invariant way, is the problem of a light quark propagation in the 
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field of the static antiquark. To simplify matter we start with the Gaussian 
correlators for gluon fields and derive the selfconsistent equations for the 
light quark propagator (with a string effectively connecting it to the static 
source). We show that CSB occurs due to the string (linear confinement), 
which shows up in the fact that effective interaction becomes Lorentz scalar. 
We also check the limit of heavy quark mass and demonstrate the usual linear 
potential in this case. As another evidence of CSB the chiral condensate is 
computed and shown to be nonzero. 

Those results were derived actually for the case of one light quark (quench- 
ing approximation or A^c ^ oo ) and when only bilocal field correlators are 
kept nonzero (Gaussian vacuum approximation). To treat the case of two 
and more flavours the corresponding term in the effective Lagrangian may 
be studied, and one is naturally led to the equation for the Green's function 
of two light quarks. The qualitative discussion of this equations leads to the 
same conclusion as in the case of one light quark-namely, the Gaussian field 
correlators which ensure linear confinement for static quarks, also yield for 
two light quarks interaction kernels growing linearly at large distances and 
produce CSB at the same time. 

Finally one could ask the question: what the effect of higher (non Gaus- 
sian) field correlators will be on the stated above results. It is argued below 
that higher correlators of even order bring about the same results as for the 
Gaussian correlators. One may then ask about the effect of the infinite sum 
of correlators, as it is in reality. The answer is that when higher order cor- 
relators only renormalize the string tension and do not make it vanish, all 
the main features of CSB hold, in the opposite case, however, additional 
investigation is needed. 

The paper is organized as follows. In the second section the general form 
of the effective Lagrangian for light quarks is given, and the equation for 
the Green's function of one light quark in the field of the static source is 
derived. In the third section this equation is studied and the conditions for 
the CSB are derived, which are satisfied in the presence of the Gaussian field 
correlators. 

The fourth section is devoted to the study of equations for the quark 
propagator in the limit of small correltion length T^, when the kernel of the 
equation becomes quasilocal. 

In the fifth section the powerful relativistic WKB methods are used to 
solve equations and calculate the kernel at large distances. It enables one 
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to calculate the chiral condensate in the limit of small Tg. This is done in 
section 6. One discovers there that parametrically the chiral condensate is 
proportional to c/T^, where a is the string tension, and therefore diverges in 
the "string limit", when a is constant and Tg vanishes. 

Therefore some additional care is needed to get the chiral condensate, and 
exact equations are written explaining these difficulties. 

The section 7 connects the chiral condensate to quark zero modes and to 
the field correlators. It is shown, that confinement occurs due to the corre- 
lator, which is expressed via monopole currents and the latter are connected 
to zero modes. 

The contribution of higher field correlators to the kernel of the equations, 
and finally to confinement and CSB is discussed in section 8. Discussion and 
prospectives are presented in section 9. 

The paper contains 5 appendices. Appendix 1 is devoted to the gauge 
-invariant derivation of the effective Lagrangian, used in the main text. In 
Appendix 2 another term in the Gaussian correlator is studied which does 
not ensure confinement, and it is shown that it also does not ensure CSB. 
Properties of the kernel, containing the confining correlator D{x), are inves- 
tigated in detail in Appendix 3. Expansion and corrections to the kernel in 
powers of Tg are given in Appendix 4. Finally in Appendix 5 the limit of 
large mass m is investigated in detail. 

2 Derivation of the effective Lagrangian for the hght 
quark 

To make discussion of confinement and chiral symmetry breaking fully gauge 
invariant, we consider the gauge invariant physical amplitude - the Green's 
function 5 of a light quark of mass m in the field of a static antiquark placed 
at the origin. The propagator for the latter can be taken as (we consider 
below Euclidean space-time). 

SqiT) = exp ig A^{r = 0,T)dT (1) 

The Green's function SqQ can be written as an integral 

(2) 
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Since S is gauge invariant, one can choose any convenient gauge for 
and our choice will be the modified Fock-Schwinger gauge, introduced in [14], 
namely: 

^4(r^=0,r) = 0, nA:(r,r)=0, i = 1,2,3. (3) 

In this gauge Gq reduces to the factor and one can now consider 
the integration over DA in (2) as a statistical averaging process and use the 
cluster expansion [15], a,b - color indices 

= exp{ j d^x'^^lx)^^ < A^ > 'ip{x)+ (4) 

1 . 

+- j dxdyij^[x)-i^ij[x)^^[y)-ir,i){y) < Ai,{x)Aj,{y) > +...} 

We have denoted the higher order cumulant contribution in (4) as . . . and 
shall disregard it for the most part of the paper, coming back to it in section 
8. The first term < > vanishes due to the gauge and Lorentz invariance 
of the vacuum, while the second can be expressed through the field strength 
correlators < Fny{u)Fx(j{u') > in the gauge (3) as follows 

A^iz, Z4) = ' duiFi4{ui, Z4) (5) 
Akiz.Zi) = r a{u)duiFik{u,Z4), a{u) = — (6) 

Zi 

One can easily see that the representation (5-6) satisfies condition (3). 
Using (5,6) one can rewrite the average < AA > in (4) as 

< Afiz)A;iiw) >= ^ l'du,a,iu) [ du',a,,{u')x 

X ^ < Fi,{u)F,,iu') > (7) 

where ab^cd are fundamental color indices and we have defined 0^4 (w) = 

1, Q;fc(w) = «(«), k = 1, 2, 3. The gauge-noni variance of the correlator 

< FF > in (7) is only apparent and one can introduce a factor, equal to 

unity in the gauge (3), which makes the correlator explicitly gauge-invariant, 

namely 

< F{u)F{u') >=< F{u, Z4)*F(w', W4)'^+ > (8) 
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where ^ is the product of 3 parallel transporters 

"if = ^{uz^, oz4)^{oz4, ow4)^{ow4; u'w^) (9) 

and 

$(P, P') = P exp ig J^^ dv^A^{v) 

where the path-ordered contour integral $ is taken along the straight line 
connecting the points P and P'. 

Since in the gauge (3) one has = = 1, we shall below omit those 
factors. For the correlator < FF > one can use the parametrization suggested 
in the second entry of [1] 

< F,^{u)Fi>^.[u') >ab= SabiS^^'6^^' - 6,^,^6^^ ^)D{u - u) + A^^) (10) 

where A*^^^ is proportional to a full derivative, its exact form is given in 
Appendix 2. 

In what follows we shall consider mostly the term D in (10) since it con- 
tributes to the string tension, while A'^^^ does not. Namely using (10) it was 
obtained in [1] that the string tension a - the coefficient in the area law of 
the Wilson loop, < W{C) >= exp^—aarea) is equal to 

a = -j_^SuD{u) (11) 

Thus a characterizes the confinement of static quarks, and our goal is to 
understand how the dynamics of light quarks is expressible through a and 
whether a correlates with CSB. 

The role of A^^^ is clarified later in section 3 when we discuss the de- 
confinement phase transition, with most details contained in the Appendix 
2. 

Keeping only the term D in (10) and neglecting higher order correlators 
like < AAA >, one obtains in (4) the following effective Lagrangian for the 
light quark: 

^effii^^i^) = j il)^[x)[—id — im)ip{x)d'^x-\- 

d'^xd'^yiip+ix)'yf,ipbix))iip^iyh„'iJaiy))x (12) 

where we have defined 

J£f {z,w) = duia^{u) du'i,a^'{u')D{u - u') (13) 
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and = 1 for ^ = 4, and = a{u) for = 1, 2, 3. 

The superscripts (^u, jj.') in J in (13) enter at zero temperature only through 
a^, and will be important for us in section 3, when we discuss the decon- 
finement phase transition. 

In what follows we disregard the perturbative contributions to >Ce//, since 
they have nothing to do with CSB. The mass m is supposed to be defined 
at the typical hadronic scale of 1 GeV, and we shall not be interested in its 
evolution to lower scales. 

From the effective Lagrangian (12) one can easily derive the equation of 
Dyson-Schwinger type for the selfenergy part, which we shall denote by M 
and the qQ Green's function S. This is done in the same way, as in the NJL 
model [16], since the structure of the Lagrangian (12) is similar to that of 
NJL however nonlocal. 

The main essential difference is the presence of the string , connecting 
the light quark to the static source, this part is concealed in J, Eq.(13) and 
therefore the selfenergy part M is actually not the set of the one-particle- 
irreducible diagrams, but rather the qQ interaction kernel. In what follows we 
shall replace in (2) the factor Gq by unity and the resulting Green's function 
will be denoted S. 

In the configuration space the equations for M and S are readily obtained 
from (12) noting that in the mean-field approximation one has to replace a 
pair of ?/^, iIj^ operators in (12) as 

Mx)i^biy) ^< Mx)i^biy) >= NcSix,y), (i4) 

and finally one obtains 

iM{z, w) = J/f (2;, w)j^Siz, w)j^dik - JltlkS{z, w)^i, (15) 

(-i4 - im)S{z, w)-i j M{z, z')Siz, w)(fz' = 6^^\ z - w) (16) 

The system of equations (15-16) defines unambiguously both the interac- 
tion kernel M and the Green's function S. One should stress at this point 
again that both S and M are not the one-particle operators but rather two- 
particle operators, with the role of the second particle played by the static 
source. It is due to this property, that S and M are gauge invariant operators, 
which is very important to take confinement into account properly. Had we 
worked with one-particle operators, as is the habit in QED and sometimes 
also in QCD, then we would immediately loose the gauge invariance and the 
string, and hence confinement. 
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3 Properties of the selfconsistent solutions 

To study equations (15-16) it is convenient to go over to the momentum 
representation. Since however the kernel J[z^ w) is not decreasing at large 
distances at |£| ~ ^ oo, we shall use at the intermediate stage the cut-off 
factor, multiplying J(z, w) by the factor exp[—a^{nP'-\-z^))^ and letting finally 
a to go to zero. Being initially a purely technical device, this trick appears to 
be of more significant contents, since it will automatically separate dynamics 
of large distances (r ~ and relatively small distances (independent of 



To make our expressions more transparent, we also assume for the corre- 
lator D{u) in (13) the Gaussian form 



We note now that due to (5), (6) the integration in (13) is done at the constant 
value of the Euclidean time component. Therefore it is convenient to treat 
the 4-th component separately, writing 



a). 



D{u) = D{0)e 




(17) 




(18) 




(19) 



and equations (15-16) have the form 




X [Ju\q.qhfiSip^,p- q,p - q)^! 



(20) 



{p4-\- p - im)S{p4,p,p} - i j 7^-^M(]?4,f,g)5(i?4, -g,^) 




= {2T:f5{p + 0) 



(21) 



where we have defined 




(22) 
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and 

K^^iz,w) = I dt I dt'a^it)a^it')e (23) 

with a^{t) = 1, ai{t) = t,i = 1,2, 3. 

In derivation of (22) we have used (11) to express D(0) through a. 
The integral over dzdw in (22) can be done easily. 

Next we introduce dimensionless momenta Q,P etc. instead of q,p as 
follows 

Q, = -^,P, = ^,P,= P±, (24) 
y/a y/a yja 

Similarly one introduces in (23) dimensionless variables r, r' instead of tt' 
connected through 

{t,t')=Tg^{T,T') (25) 

Now in the kernel Jik the a dependence can be explicitly written as 

Jt^ = ''-^ffi(Q<Q') (26) 

Jik = -fnkiQ.Q') (27) 

here fik are dimensionless functions of order one, when their dimensionless 
arguments are also of the order of one, fast decreasing at infinity of Q,Q' 
and finite at small Q,Q'. The important fact about fik is that they do not 
depend on a in the limit Tgy/a — ^ 0, which is solely of physical interest. 
For small Q, Q' one has 

mO,0) = lSi^^'" (28) 

From the dimensional analysis of equations (20), (21) one can derive the "a 
- dimensionality" of M and S, namely: 

M{j)^,p,p') = a-^/^M{P^, P, P') (29) 

S{p^,p,p) = a-^S{P4,P,P') (30) 

Insertion of (29), (30), (26), (27) into (20), (21) yields dimensionless equations 
for M and 5 

(P474 + Pm - ^)S{P,, P, P')- (31) 



-./||^^^5(A.-Q.P') = (2.f.'^'(P + P') 



X745(Pi, P-Q.P'- Q')74 + QOT^'^CPi, P-Q.P'- Q'hi 

- aT^fSiQ, Q')lkS{Pl P-Q,P'- Q'H} (32) 

In the rest of this section we shaU analyze the system of equations (31), (32) 
in the physical limit, when aT^ 0, while aTg is kept fixed. 

To understand better properties of solutions of (31), (32) it is convenient 
to consider first the limit of large mass m. In this case one expects to get the 
usual Lippman - Schwinger equation corresponding to the linear potential 
between two heavy quarks. 

The limit of the large mass m has its own peculiarities which can be easier 
seen in the configuration space. Hence we come back to equation (15), where 
one should insert the large-mass propagator Sm{z, w) 




27r (27r)6 pi + 



(33) 

obtained in the lowest order in (^) from (21). Thus one has the usual form 
for /i4 > 

SUhi^z.w) = '-^^^^0{lH)d^''\z- w) (34) 

and we have defined = z^ — w^ and have factored out the term exp{—mh4), 
which contributes to the total factor exp{—mT) of the overall qQ Green's 
function. 

Now we introduce Sm or the r.h.s. of (15) and compute the kernel Jik 
given in (13). One has 

JS'ihA, £, z) = D{0)e~^ I' I' a^{t)a^it')dtdt'e~^^'~''^\iZk (35) 
At large \z\ one obtains 

Jf/(/i4, z, z) = D{())e-^hTl^^^z,z, 1^ I 1, 2, 3 ) ^^^^ 

Expressing now i^(0) through the string tension a via (11) one obtains for 
M{z,w)] 

1 

M(/i4,i',^) =e"^i:'(O)0FT5|f|(2 + ^^^)(5^3)(£-w;)^(/i4) (37) 
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The equivalent static potential U (z, w) obtains when one integrates over the 
relative time in the limit Tg 0, and keeps only the (+, +) component in 
Lorentz indices (or equivalently puts 74 = 1) 

U++{z,w) = J dh4M{h4,z,w) = a\z\S^^\z- w) (38) 

Thus one indeed obtains the standard linear potential for the heavy quark, 
which is Lorentz scalar. 

We now turn back to the equations (31), (32) and look for the hmit a ^ 0. 
To this end we remark in (32) that while P, P'Q, Q' are confined to the finite 
limits independent on P4 and P4 are not and can be of the order of 
therefore one can neglect Piji on the l.h.s. of (31), thus writing it as 

p _|_ i{m+M) 
^4 H a 

Inserting this in (32) and integrating over (iP4 in the limit aTg one 
obtains 



iM{P,, P, P') = '^l ^^(2^fi^'\P + P'-Q-Q') 



(40) 



2 -I (27r)< 

[(m + M)^]-i/^(m + M)/«(Q,Q') 

Here M on the r.h.s. is actually an operator and its dependence on momenta 
in the integrand should be properly written, as well as the definition of the 
power (—1/2). We would like to make several comments on Eq.(40). 

i) First of all one can see that the r.h.s. of (40) does not depend on a, so 
that all expressions can be considered as in the physical limit. To see it more 
clearly, one can use (30) to express the physical Green's function S through 
the mass M as follows 

S{Pi,P,p') = ^^,{2^)V'>(p + if) (41) 

pi + 

and one can see that M enters the physical propagator, and is the real physical 
quantity. 

ii) Secondly, only the A4 component (which generates f^^) contributes 
in the limit aT^ 0, while Ai,i = 1,2,3 give vanishing contribution. In 
terms of field strengths it means that only color-electric correlators < EiEj > 
contribute to M in the physical limit a — ^ 0, while color-magnetic do not. 
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This fact seems to be connected to our choice of gauge, since in this gauge 
q;^=4 = 1, while a^=i = t and the latter yield suppression of the contribution 
of Ai^ at large distances. However, it is clear that the gauge corresponds to 
the physical situation, when the static quark has the world-line along the 4-th 
coordinate, and therefore the string which is formed between the light quark 
and the static one, is evolving alone the 4-th coordinate, hence it is done by 
the electric field. In the last section we consider this point more closely, since 
it explains why CSB disappears together with electric correlators D(x)(10) 
at the deconfinement phase transition. 

iii) The most important feature of (40) is that it exhibits a finite scalar 
solution for M even in the limit of vanishing quark mass m. This solution 
exists for any finite value of crT^, and is the consequence of behaviour at large 
distances (corresponding to the limit a ^ 0). 

For less singular behaviour of the interaction kernel Jik at large distances, 
such as (27) instead of (26), the selfconsistent solution (40) does not exist. 

iv) In (40) we have retained only terms of the correlator D, neglecting Di. 
The role of Di is discussed in Appendix 2, where we show that it does not 
give selfconsistent solutions in contrast to the case of D. 

4 The local limit of selfonsistent solutions. Confine- 
ment and CSB for light quarks 

Our bassic equations (15), (16) are nonlocal in time because of the integral 
over dz'^ in (16). This nonlocality and the parameter which it governs can be 
handled most easily, when one uses instead of M(^z^z')^ S{z,z') the Fourier 
transforms. 

Siz, - zl z, ^) = I e*^^(^-^^)5(p4, (42) 

and the same for M{z, z'). Then instead of (15), (16) one obtains a system of 
equations 

(p4 — idz — im)S{p4, z,w) — i J M(p4, z, ^)S{p4, z', w)d^ = 6^^\z — w) (43) 
X [ Jf/(z, wh.Sip',, z, wh,6ik - JihkSjil (44) 
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where 5^ = 7^^, J is defined in (13) and we have factored out the time- 
dependent exponent of D{u), using the representation (17). (For any form of 
correlator D{u) the main result below, Eq. () remains true, but the correc- 
tions to it are dependent on the shape of D{u) and are displayed in Appendix 
4). 

All dependence of M on as can be seen in (44) is due to the factor 
exp[—{p4 — p'^yTg] and disappears in the limit when Tg goes to zero, while 
the string tension a ~ D{0)T^ is kept fixed. This limit can be called the string 
limit of QCD, and we shall study its consequences for equations (43), (44) in 
this section, while in the appendix 4 corrections to this limit are considered. 

So in the string limit, with M independent of ^34, let us consider the her- 
mitian Hamiltonian 

H^n = (-^ + pm)M^) +P f Mip, = 0, z, ^)M^)d'z' = e^M^) (45) 

where 

with eigenf unctions V'n satisfying usual orthonormality condition 

The Green's function S can be expressed as 

5(p,f,,T=S«S«M (46) 

n P474 - ^£n74 

Inserting (46) into (44) one is met with integrals of the type: 



/- 



00 dpA e^^'P'^^^'^^^^^a i 

Note, however, that the result depends on the boundary conditions. If, e.g., 
one imposes the causality-type boundary condition, then one obtains 

We are thus led to the following expression for M in the string limit 

M(p4 = 0, £, w) = y/7rTg[Jli''{z, w)^^K{z, w)^4l^ - J'i{z, w)^kMz, W)j4li\ 

(48) 
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where the definition is used 



k{z, w) = Y.^n{z)sign{£n)il)niw) (49) 

n 

Let us disregard for the moment the possible appearance in M of the 
vector component (proportional to 7^, // = 1, 2, 3, 4) and concentrate on the 
scalar contribution only, since that is responsible for CSB. Then one can look 
for solutions of the Dirac equation (45) in the following form [17] 




Gn{r)QjlM 

iFn{r)^ji'M 



(50) 



where V = 2j — I, and introducing the parameter K,{j, I) = {j + ^)sign{j — I), 
and replacing M by a local operator (the generalization to the nonlocal case 
is straightforward but cumbersome, the final result (53) is not changed in the 
nonlocal case). We obtain a system of equations [17] 

^ + ^Gn-i6n + m + Mir))Fn = 

+ f i^n - (£n - m - M{r))Gn = ^^'^ 

Eq.(51) possesses a symmetry {sn^Gn, F^, k) ^ (— F„, G^^, — /t) which 
means that for any solution of the form (50) corresponding to the eigenvalue 
there is another solution of the form 




Fn{r)^jl'M 

iGnir)^jiM 



(52) 



corresponding to the eigenvalue (— £n)- 

Therefore the difference, which enters (49) can computed in terms of Fn, Gn 
as follows 

A(^,y) = — E [Gn{x)Gi{y)nM^WjiMiyy 

•^y njlM 

-Fn{x)F:{y)^,i,M{mji'Mm + — E {Fn{x)Gl{y)x 

•^y njlM 

X nj,M{x)n'jiM{y) - Gn{x)F:{y)njiUx)n*ji'Mm (53) 

The expression (53) serves to display the appearence of the 74 factor in the 
main term, which as one can undestand from (48), produces the scalar con- 
tribution for M(z, w). For quarks of heavy mass the sum (53) reduces to the 
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^-function term, which can be most easily seen in the simphfied example, 
taking M to be constant. In this case one can write 

One can easily see from (54), that for (m + M) oo one has 

limAm+M^oo (^, y) = 74(5^^^ {x - y) (55) 
We now turn to the Jikiz^ £), (13), which can be rewritten as 

Jt^[z, z) = ZiZu I' I' dtdt'D{z{t - t')), 



Jll{z, z) = ZiZk I' tdtt'dt'D{z{t -t')). (56) 
For the Gaussian form (17) one obtains at large \z\ 

4Mil ^ ^) = ^2T,^D(0), = \^2T,ViD(Q) (57) 
Inserting this into (48) one gets for M{z, z) at large distances \z\ 

M{p^ = 0, £, w) = a\z\{2 - ^-^)6^^Kz- w) (58) 

Thus one obtains seemingly a local Dirac equation for and S with the 
mostly scalar kernel, containing linear confinement. 

However, the kernel (58) does not correspond to the Dirac equation for a 
quark of large mass (m + M) (the latter property was used to replace A{x, y) 
by the (5-function, as in (55)). 

To see this, one should keep in mind that the kernel (58) is the limit 
P^Tg ^ of the full kernel M(p4.z. ^Z;), Eq.(44). while for the large mass 
m the effective j94 is arround m and therefore also large and the local Dirac 
equation obtains in the hmit, when one uses instead of equations (43), (44), 
the set (15), (16), where in the lowest approximation in M/m one inserts in 
(15) the free propagator for the heavy particle m, Eq. (34). Details of this 
derivation and results are given in Appendix 5. 

One can see that the resulting Dirac equation has the effective kernel (38), 
with the expected coefficient (j\z\, in contrast to (58). Therefore we shall see 
that the expression (58) is valid only for light quarks, mTg <C 1, and for the 
case when A(i', w) reduces to the 5 - function term. To see when it is possible 
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one can use the quasiclassical approximation to calculate A(z, w) and we do 
it in the next chapter, while we conclude this chapter with the analysis of 
the qQ spectrum, resulting from the spectrum of the Hamiltonian (45). 

There is an important point, which one should have in mind concerning 
the Dirac equation. As we have discussed above in this section, the Dirac 
equation with the scalar interaction produces the spectrum symmetric under 
reflection Sn — and it is clear that negative energy states are not present 
in the spectrum of the heavy-light mesons in reality. To resolve this paradox, 
let us come back to our Green's function 

SqQ = (rS'(x4 - 1/4, z, y)rSQ{x4 - Va)) (59) 

The masses of the {qQ) system are obtained from the asymptotics of SqQ at 
T = X4 — y4 ^ oo, i.e. SqQ ~ exp{—EnT) one should put in (59) Sq as 

SQih) = ^[l^e(h)e-"^'' + L^e{-h)e"«h] (60) 
and considering positive T = X4 — y^, one can express GqQ as 

n ^ ^ J P4-i£n 27r 2 ^ 

= E^^(r74|n > e{£n) < nle-^^Q+l""!)'^ ^^"^'^ ) (61) 

n 2 

Hence one can see that only positive values of contribute to the mass of 
the qQ system, namely 

MqQ = MQ+\enl 

Using (48), (49) and (45), one can write a nonlinear equation for eigenfunc- 
tions ipn, namely 

d r ~ 

C^i '' k 

X747m - <4^(^,w;)7jEV^/fc(i)si5rn(£fc)?/;J(w;)747i]V^n(w;)rf"^w; = 



A; 



where we have defined 
and J-^ is given in (13). 



= enipniz) (62) 
Jj^"^ = V^T,4t^ (63) 
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5 Quasiclassical solution of the selfconsistent equation 



In this chapter the quasiclassical analysis of equations (62) will be given. The 
nonlinear part of the kernel M, Eq. (48), is A(x, y) which can be expressed 
through solutions ipn as in (53). Our primary task is now to calculate A(x, y) 
in some reasonable approximation, having in mind to improve it at later step. 

Therefore we start calculating A(x,y) in the form (53) using for 
quasiclassical solutions of the local Dirac equation with the scalar potential 
U = err. As one will see the resulting A(x, y) is the quasilocal object, tending 
to the 6{x — y) at large \x\. Being inserted in (62) or (48), A{x,y) would 
indeed generate a local kernel M(z, w), producing the local (at large r) Dirac 
equation with U (r) tending to err at large distances. Thus one justifies a 
posteriori the initial choice of Gn, at least for large r. 

Coming back to calculating A(x, y) for U = err, we follow the method of 
[18] and write the Dirac equation as follows 

drr^=n^^-^^=[F)^^=[m'iu-eU ^ ) ^''^ 

We have kept the vector potential V in (64) to make our consideration as 
general as possible. 

To do the quasiclassical expansion, one writes [18] 

00 



Ip = ipexp j y{r)dr, p = ^ h^(p^^ 



n=0 



y{r) = ^y-i{r) + yo{r) + %i(r) + ... (65) 
and obtains the system of equations 

{D - ?/_i)^(") = 

(D-^_O^W = |-^W + yo^W, 

(D - ?/_i)^("+i) = -f + E ynW^ (66) 

dr k=o 

In what follows we keep in the lowest approximation y-i{r),yo{r) and ip^^\ 
and obtain 



y-i = Xi = ±^ 



k2 



^ (m + t/)2 + — - (£ - Vy = ±q (67) 
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A' 



+ 



U'-V 



K 



{U' - V] 



2Xi 2{m + U + £-V) 2X,r {m + U + e - V) 



^ = A 



m + U + e - V 



= A' 



Xi — - 

m-\-U - e-\-V 



, (68) 



(69) 



and 



ifj^^ = iff^exp f {y^i\{r')dr' + yf{r')dr') (70) 

where the index z = +, — refers to two possible solutions of the matrix equa- 
tion {D - y-i)(f^^^ = 0. 

From (67) one can find three different regions on the line < r < oo ( for 
U > \V\). Taking for simplicity m = 0, V = 0, one has from ^ = 0, 



2^^ 



(71) 



Thus for g2 > 2aK ( and this holds for all levels, see below) one obtains two 
turning points r± in (71). Between these points Xi is imaginary and this is 
the classically allowed region with the momentum p[r) 



p{r) = 



(g - "\/)2 - — - (m + uy 



(72) 



The quasiclassical solution in this region can be written in analogy with the 
corresponding solution in [18, first entry Eq. (1.12)] 



G = C 



m + U + £ - V 
p{r) 



sinOi; F = CisgriK 



e — V — m — U 



where we used the notation similar to that of [18] 



p{r) 



sin62 (73) 



e,{r) = l_(p + —)dr' + -, 0^(r) = Jjp + —)dr' + - 



w[r) = —- 



U'-V 



2r 2m + U + e -V 



(74) 

(75) 
(76) 

2r 2m + t/ + g-y ^ ^ 

Here prime denotes the derivative in r; putting U = one recovers 
formulas(1.12-1.14) of [18]. 



w[r) = 



U'-V 
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As a next step we turn to the quasiclassical determination of energy eigen- 
walues Cn- Using the Bohr-Sommerfeld equation 

pp{r)dr = 7r{n + ^), n = 0,l,2..., (77) 

one obtains for m = 0, y = 0, with p(r) from (72) 

el = Mn + ^^) = 2a(2n + j + ^) (78) 

Comparison of the numbers obtained from (78) with exact calculation of 
Dirac equation in [17] reveals that (78) is a very poor approximation for 
the real spectrum, and, moreover, it is qualitatively incorrect, since Sn in 
(78) does not depend on the sign of k, i.e. on the spin-orbit interaction 
of the light quark. Happily this problem was already treated in [18] for 
Coulomb interaction, and the authors have proved that it is legitimate to take 
into account spin-orbit interaction quasiclassically; they suggested another 
eigenvalue condition instead of the Bohr-Sommerfeld, which can be deduced 
from the solutions (73), namely 

rip + —)dr = 7r(n + ^), n = 0, 1, 2, ... (79) 

Jr- pr 2 

where w is defined in (75). Doing the integrals approximately in (79), one 
obtains the following spectrum instead of (78). 

9 . / 11 1 + sqriK, K,a ,^ , ei ^f,K.o-.r,\ , 
el = Mn + 2 + 2^ + " 2^(0-38 + ^n-^)) + ((^)^) (80) 

One can visualize the appearence of the last two terms, depending on the 
sign of which produce the spin-orbit splitting of the levels. The form (80) 
is exact for small (^)^, and for lowest levels the accuracy can be tested by 
comparison with the exact solution of Dirac equation [17]; two sets of 6 lowest 
levels coincide within 5%, whereas for (78) the accuracy could be as bad as 
25 -i- 30%. Note however, that for high excited levels, where n and /or \k\ 
are large, the correction terms in (80) are constant, and therefore relatively 
unimportant as compared with growing n and j. Therefore in what follows 
in summation over high excited states we shall use the simplified form (78) 
instead of the corrected one, Eq.(80). 

We turn now to the quasiclassical calculation of A(x, ^), Eq.(53), and to 
this end we represent it in the following form 

A(f , y) = /3Ai(f , y) + i75A2(f , y) (81) 
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where, e.g., Ai is a 2x2 matrix, /i, /i' = ±| 



•^y njlM 

- Fr,{x)F:{y)^^^,^{x)^%{y)}, (82) 

and the sum is over positive energy states only; Gn^ depend not only on 
n, but also on j, / but independent on M. 

We shall be interested mostly in the large values of n, / in the sum (82), 
since those terms will form a tempered (5-function of the type 5^^^ {x — y). 
Hence one can neglect in 6*1,6*2, Eq. (74) the constant terms ^ as com- 
pared to the growing contribution / pdr' . Now in p{r), (72), one can replace 

by at large /, and as a result both Gn and F^ in (82) do not depend on 
j (with / fixed) and one can sum up over j, M 

T.^jiMi^)^jiMiy) = (^wE^/,M-M(^)>l>-^(y) = 

jM M 



Hence A'/'^ is diagonal matrix A^'^ = ^uu'Ai with Ai equal to 



(83) 



2/ -u 1 

Aiix,y)=Y.G!inJ)^—Pi{cose,y){[Niix)Niiy)-N2ix)N2iy)]sme{x)si 
n,l 47r 

where we have defined 



(84) 



Ni(x) = 



1 



£n + m + U{x) — V{x) 
p{x) 



■,N2{x) = 



£n — m — U{x) — V{x) 
p{x) 



(85) 



and 6{x) = /;?^ pdr. 

The normalization constant Ci(n, I) entering (84), can be estimated as 

C( dr{ ^j—^ sin^Oi H 3111^92} = 1 (86) 



p{r) 



p{r) 



For large /, n the functions sin6i{r)^ sin62{r) are fast oscillating and hence 
one can replace 

(87) 
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and one obtains foi V = 0, 

C- = ef^4- = f^; C!in,l) = ^ (88) 

Jr- plj-j 2(7 er^TT 

Consider now in Ai(x, y) both x, y in the region 

r_<Ca:,y<r+, — y|<Ca:,y (89) 

Then one has 

N,ix)N,iy) - N2ix)N2iy) ^ ^^"^^^^""^^ (90) 

p[x) 

and ^ 

sin9{x)sin0{y) = -[cos{6{x) — 9{y)) — cos{9{x) + 9{y))] ^ 



1 

-cos 



/ \ £^ — cr^r^ T^dr + ... (91) 

where dots stand for a fast oscillating term, which will not contribute to the 
final answer. Hence one obtains for Ai, 

Ai f , y = — E -^Piicose^y) ^^cos / pir)dr (92) 

One can compute the integral under the cosine 

1= p[r)dr = -{p{x)x-p[y)y)- — {a{x)-a{y))-—{(3{x)-(3{y)) (93) 
Jy 1 4(7 2 

where p{x) is given in (72) with 1/ = 0, and 



^2 _ 2^2^2 e"^ -H- 

a{x) = arc sin ^====, /3 (a;) = arc sin f = (94) 

In what follows we replace the sum over n in (92) by integration over £, with 
/ (or fixed and neglect in (93) terms containing which are small at 
large n. 

We also introduce instead of e the variable r as follows 

e = (JXt; t > l;p{x) = crxy/r'^ — 1 (95) 
From (80) one has asymptotically 
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Expanding a{x,y) in (93) around x = y one has 

I{x, y) = p{x){x — y) + 0{{x — y)^) = ax{x — y)y/r'^ — 1 (97) 
Finally one obtains in (92) 

X I 'K All V'T — 1 

= — E ^-^Pi{cose,y)K,{a) (98) 

where 

a = Gx\x — y\ (99) 

and Kq is the McDonald function. 

Now one can perform the sum over I in (98) to get 

E(2/ + l)Pi{cose) = 4(5(1 - cose) (100) 

Thus one obtains for large x, y and |a; — y| <^ x^y the result, which we write 
in the form symmetric in y 

Mix, y):,^y = ^ K^{a^fxy\x - y\)b{\ - cosd^y) (101) 



One can test that Ai(x, y) is a smeared b - function, normalized to 1, with 
the smearing radius 

ro ~ [a^JxyY'' ~ (crx)""^ <C mm(^^,x) (102) 

The normalization can be checked doing the integral: 

/ Sykx{x,y) = [ dcose6{l - cose)dipy^^Ko{a) (103) 

writing 

, a , , da ^ .X 

?/ = a;± — , dy = ± — 104) 
ax ax 

and taking into account that 

J d cose6{l - cose) = ^ (105) 

TOO VF 

daKo{a) = - (106) 
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one obtains 

Jd'yAi{x,y) = l (107) 

Thus Ai{x, y) is & normalized smeared ^-function, which is "focusing" the 
nonlocal interaction kernel M(x,y), Eq. (48), at least at large distances x,y 
into a quasilocal kernel, linearly growing at large x ^ y. 

At the same time another part of A{x, y), namely A2{x, y), (see Eqs. (53), 
(81)) can be written as 

Afix,y) = - E {Fnix)G:iy)n^,MimjMiy)- 

njlM 

- G„(x)K(^)n^,^(x)Q*i;;^(^-')} (108) 

Doing the same procedure with A2 as was done with Ai, and to this end using 
the relation 

^ji'Mix) = -(o^nx)Mi^O^zM(^) (109) 

one obtains 

-1 



Af [x,y) = —XCt{n,l){N2[x)Ni[y)sme2[x)smei[y)^^^^^ 



2/ + 1 

-Ni{x)N2{y)smei{x)sme2{y){afiy)^^>}— Pi{cose) (110) 

One can see from (110) that A2(x, y) is odd with respect to exchange x ^ y 
and therefore vanishes at x = ^. Hence A2{x, y) has no local limit and gives 
no contribution to the long distance linear interaction A{x^x). 

We shall not be interested in A2 for the rest of the paper, however it may 
contribute to the finite range nonlocal part of the resulting kernel M(£, y). 



6 Calculation of the chiral condensate 

Our method allows to calculate the chiral condensate - the characteristics, 
which does not depend on the presence of the static quark. To this end we 
consider the euclidean condensate, which due to the definition (14) is 

< il)^{x)il){x) >= NctrS{x, x) (111) 

At this point it is important that we take in S{x, y), x = y = 0, thus putting 
the initial and final point on the trajectory of the static quark, and not at 
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the point x = y far from origin. The reason is that with our choice of the 
gauge, (5), (6) and in the Gaussian approximation (12) the string from the 
hght quark extends from the point x, y necesarily to the x = and therefore 
the world sheet of the string is not the minimal surface, but rather a cone 
with the vertex at the point x = 0. The choice S{Ox = 0, ^ = 0) makes 
the area inside the characteristic trajectory of S as minimal as possible, thus 
greately diminishing the contribution of non-Gaussian correlators. 

The latter collectively are responsible for giving the minimal area contri- 
bution for any choice of the point zq in the Fock-Schwinger gauge condition. 

To proceed one can use our definitions (46), (47), (49) and (81) to write 

< V^+(0)V^(0) >= trSih^ = 0, 0, 0) = 

/t^^r5(p4,0,0)=iAr(0,0) (112) 
Using (82), one can rewrite (112) as 

The analysis of the behaviour of Gn{x), Fn{x) as solutions of the Dirac equa- 
tion with the linear potential at small x was performed in [17], where the 
following properties were found 

G = Ax^''\ + ...,F = Bx^''^, (114) 

Thus in (113) the contribution comes only from the states \k,\ = 1, i.e. j = 1/2 
and I or I' = 0. 
In the first case 

l)/t = -1, / = 0; A ^ 0, B = (115) 

In the second case 

2)k = 1J = 1;1' = 0;A = 0,B^O (116) 

Denoting 

A{n, K = -l) = B{n, k = +1) = B+ (117) 

and heaving in mind the relation between Minkowskian and Euclidean quark 
condensate 

>M=i <iJiJ^ >E (118) 
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one obtains the general expression 

Af 00 

< # >M= E [(^n)' - (Kf] (119) 

Our next step is the calculation of A^^B^ using the quasiclassical method. 
The standard matching condition connects the wave function below and above 
the turning point r = r^: ^ 



where si is finite for r ^ 0. 

Similarly for 62{r) one obtains 



From (120), (122) and (123) the the A' and B+ are 



where the following notations are used 



e+ = £(n, = +1), = £{n, k = —1) 
= Ci(n, K = ±1) 



(120) 



Denoting 

— rr Kill — rr Kill 

Oi{r) = / {p+—)dr',e2{r) = / {p + —)dr' (121) 
one has the following behaviour of 9i{r) at small r: 

- 1 7" 

9i{r) = i(|«:| + -signtijln- h ^si (122) 



1 r 

92ir) = ^(k| + -signK)ln- h ^S2 (123) 

2 2r_ 



= 5+ = ^i^e^^ (124) 



(125) 



To the lowest order in 1/n expansion one has 

II K t: tt ^/ 1 X 

^^ = i"i-2i^2 = ' + 4i^+°^;?^ 

S2 = \A - = 1 - ^ + O(^) (126) 



^ Recently an important correction was found [19] to the standard nonrrelativistic matching condition, 
which is not taken into account below, since at large n it is small. 
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Since for large n one has 



(4)2 = Mn + l±J + 0(i)) (127) 



and due to (88) C± = (^)^/^ and r_ = ^(l + 0(^)), the difference {A-f- 
(B+)^ vanishes at large n in the leading order, namely one obtains 



TTV^ + 1 

In (128) the factor in brackets occurs from subleading terms in (the term 
(-1)) and from the corrections si, S2 (the term tt). However the corrections of 
the order of 1/n in are not taken into account in (128), moreover, the next 
order quasiclassical expansion terms, i.e. yi{r) and (f'^'^^ may contribute to 
the 1/n correction to the coefficients A~, B^, and the problem of computing 
the difference (128) becomes rather complicated. 

With all corrections included, the difference {A~y — [B^Y is of the 
order 0(-^) ( modulo unexpected cancellations), and the series (119) di- 
verges. Therefore we must try to approach the problem of computation of 
5(0, 0) ~< V^V^ > from another side. To this end we first consider the prob- 
lem of computing V^ri(O) for the nonrelativistic Schroedinger equation with 
the local potential lJ{r). Following [20], one has 

i,„(o)i^ = i < nr) >- i^-^^^^ im 

In particular, for the linear potential U{r) = or, |?/'„(0)p = ^ and does not 
depend on n. For the quadratic potential one has |'0„(O)p ~< r >nn'^ \fn. 

As one can see in (128), this case is similar to ours, i.e. the effec- 
tive Schroedinger potential for the Dirac equation with linear interaction 
is quadratic. This fact is well known in the quasiclassical approximation 
[17,18], indeed one can write the equivalent Schroedinger potential with the 
energy and effective potential IJ of the following form [17] 

E = \e\ U = \U'^^. p\t) = 2{E-U), 

(P_P. + U)'il;{r) = EiPir),/! = 1 (130) 
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Taking U = ar, one immediately obtains the leading term for {A~y, (B^)'^ for 
large n, proportional to a^^'^y/n. Thus far the effective nonrelativistic theory 
with a local potential confirm our result (128) and produces the diverging 
sum as in (119). 

Let us now consider the nonlocal effective potential, U (r, r'). The straight- 
forward calculation similar to (129) yields in this case 

\M0)\' = ^ j^dr dT'y:{T){^U{Ty))yn{r') (131) 

where yn{f) is the radial part of ipnir)^ normalized as y'^{T)dr = 1. To 
understand to role of the nonlocality, let us represent U as follows 

U(r^ r') = (v(r) + v(r') f ^ (132) 

2a^/7^ 

so that in the limit a — ^ one has 

d(r — r)U(r,r) = v(r) (133) 

-oo 



For the case when v{r) corresponds to the equivalent efTective potential, i.e. 

1 

2< 



v{r) = \a'^r'^^ii = 1, see (130), we obtain 



/-oo /• _ , ^ , , r (r-r' f 



(J /"OO r I \T T ) 

= -i J drdr'y:ir)^e-^ynir') (134) 
For the chosen above potential v{r) the wave function yn{r) is 



C r 
Vnir) = -^^cos[{2n + 1)(- 



r r TT 

1 — ( — )^ + arcsm — ) — — ] (135) 



lp{r) r+\ r+ r+ 4 

where = i(n +1/2),C^^ ^, and 



\\^n~ <7^^^ 2 ~ V^^ ~ ^^^^ (136) 



For a <C r+ and large n one can approximate the integral (135) as follows 
1/ /nM2 ^2^^^ f d{^)r fd{r-r') An , , {r-r'f 

(137) 
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where (f{r,r') = 1 — ^^"'"g^^"'"'^ ^ + 0((^)^). The integral over d{r — r') yields 
the cut-off factor 



1 foo 4n (r — r'Y _4nf^^ 
—j= I dir — r')cosi — (r — r)(p)exp( ^ ) = e (138) 

and one has 



\M0)\' = —eM-^-^^') (139) 

where we have defined the effective value of of the order of one. 

Let us compare (139) with the square of A~,B^ from (124). The first 
factor in (139) coincides with that of (124) up to a factor (^), while the 
cut-off exponential factor is new and originates from the nonlocality of the 
interaction (132). It essentially cuts off the sum of (119) at the n > Umax, 

11 . X 

rimax - ^ o (140) 

We turn now back to the calculation of the chiral condensate, Eq. (119), 
with the A\ — given by (128). The last expression is obtained without 
corrections 0(l/n) to the normalization constants C^, which may change the 
numerical coefficient in (128), but cannot change the n~^/^ behaviour, leading 
to the divergence of the sum (119). 

We have noticed above, Eq. (139), (140), that nonlocahty of the effective 
potential, i.e. M(p4, z,w) causes the cut-off of the sum (119) in n, and now 
one must look more closely at the origin of this nonlocality. 

To begin with we should remember, that the expression (82) for Ai{x,y) 
and as a result the sums in (113) and (119) are obtained in the limit Tg 0, 
when according to (44) Ai and M(j94, w) do not depend on p^, so that we 
have put j94 = in M{p4, w). Let us now take the finite Tg into account and 
compute Ai(x, y) anew. To this end one can use the results of the Appendix 
4, and for simplicity we again put there P4 = 0. In this case the corrected 
A{x, y) is 

Hx,y) = EM^>^9nenfi\en\Tg)i/jiiy) (141) 

n 

where from (A. 4. 7) the function / is (the exact form is given in (A. 4. 5) 



2|£n|T, 



\en\Ta <C 1 



fi\en\Tg) = \ ' irif ; (142) 
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On can check that the presence of this function in the sum (92) and in the 
integral (98) does not change the main result: the appearance of the smeared- 
off ^-function with the same range parameter (99), but instead of the function 
^Ko(a) one obtains another smeared ^-function K(a), finite at a = and 
exponentially decaying at large a. Therefore again one obtains local at large 
r interaction - linear confinement for light quarks. 

More important is the change in the sum (119) which now should be 
replaced by 

/V oo 

< >M= E KA- f - {B:m\e„\T,) (143) 

Using the asymptotics (142) and (128), one can see that the sum is now 
diverging only logarithmically, and this weak divergence may be cured by the 
nonlocality caused by the rest dependence of M{p4, z, w). 
As a result one obtains an estimate 

< >^ _iv^^/„^ (144) 

Tg So 

where A is the numerical factor, Smax ~ l/^^- For Tg ~ 0.2fm and numerical 
factors from (128) one obtains a correct order of magnitude for the chiral 
condensate. A more detailed calculation of < > in the quasiclassical 
method is now in progress and will be published elsewhere. 

Note, that the Nc dependence in (119) and (143) is reproduced correctly 

- the chiral condensate should be proportional to Nc in the limit Nc ^ oo. 

Since a ~ D{0)Tg ~< F^(0) > Tg, one obtains for the chiral condensate 
parametrically (up to a numerical constant) 

< >= < F\0) > Tg (145) 

TT 

Thus the chiral condensate diverges in the "string limit of QCD", i.e. when 
a = const, and Tg 0. In this limit the width of the QCD string tends to 
zero [21]. In the realistic case, i.e. for finite Tg, and finite gluonic condensate, 
the < i/ji/j > is also finite, and inserting in (145) the standard value [22] 
f < F\0) >= 0M2GeV\ and Tg = IGeV"^ [23], one obtains < ^j^j 

— (250 MeV)^ const, with const of the order of one. 

To understand better the source of divergence of the chiral condensate at 
Tg ^ and to prove the finiteness of < i/ji/j > for the exact equations (15,16), 
let us look for the solution of (16) in the form 

S {x, y) = -idxfix, y) + g{x, y) (146) 
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Insertion of (146) into (16) yields for m = 0. 

- dlf{x, y) - id^gix, y) - i j M{x, z){-idj{z, y) + g{z, y))(fz = 6^^\x - y) 

(147) 

where 

iM{x, z) = J^{x, z){-id4-f4f + g + idif)+ 
+J^{x, 2;)[+3i54/74 + Sfif + id^f - 2inzdfhx - [g + idf)hzhx] = 

= iM\^dakf{x, z) + M^''\x, z)g{x, z) (148) 
To satisfy (147), f{x,y) should be singular and one can represent / as 

and 

As a result one obtains two equations 

1 '--rW, 



-dtf{x, y) = -tr J [Ml^>-iAf[x, z)^k'dk'f{z, y)+ 

+ M^'^\x, z)g{x, z)g{z, y)]d'^z (151) 

-id^g{x, y) = i^tn^ J [M^^^difix, z)g{z, y)-M^'^\x, z)g{x, z)^,ydj{z, y)]d'^z 

(152) 

One can identify on the r.h.s. of (151) the most singular term, which yields 
for / at small x, y 

f{x,y) ^ ln\x -ylix"^ + xy) (153) 
where we have the property of J^, at small x, y 

J^, - {x^y^)const (154) 

From (152) one can conclude that g{x,y) is nonsingular at small x,y: 

9[x, y) = go + giix^ + y^) + g2xy + ... (155) 

Hence 

< ijjijj >= m,g{0, 0) = m,go (156) 
is nonsingular and finite for the solution of the full equations (147). 
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In the small Tg limit, considered in the last two sections, when the kernel M 
is becoming local and yields the linear scalar potential at large distances the 
situation changes. One can do the same analysis as done above, Eqs. (146- 
152), and obtain for the local time-independent kernel M = cr |x| (5*^'^'' (z — y)) 

hm a(x,0) = -^, (157) 
NHo \x\ 

and when the time-nonlocality term is used in M, M ~ exyi^— ^^^^^ ), g{x, 0) 
is less singular g{x^ 0) ~ /n|x|, x ^ but still diverges logarithmically. These 
properties explain the divergence of the sum (119) in the approximation (128), 
corresponding to (157), and the logarithmic divergence of the sum (143) corre- 
sponding to the time-nonlocal case (with the account of the j94-dependence). 
Our conclusion is that the accurate computation of the chiral condensate 
requires exact solution of the equations (15), (16). 



7 Chiral symmetry breaking: zero modes vs field cor- 
relators 

We have seen in previous sections that the nonhnear equations (15), (16) give 
rize to the phenomenon of CSB, which reveals itself in our problem in two 
ways: i) it provides scalar confining interaction for the light quark ii) there 
appears a standard chiral condensate < ?/^?/^ >. 

A natural question arises at this point: a folklore understanding of CSB 
is that it is due to quasizero global quark modes in the vacuum. An exact 
relation [24] exists, which connects chiral condensate to the density z/(A) of 
quasizero modes in the vacuum at A < m, and in the chiral limit (m 0) 
one has [24] 

< >= (158) 

Here A is an eigenvalue of the 4d Euclidean equation for the quark in the 
vacuum field Ay^ 

- ibll)n{x) = Xni^n{x) (159) 

The density v{\)d\ is the averaged over all fields {A^} number of the states 
A^ per interval d\. 

It is a popular belief that the quasizero modes necessary for CSB due to 
(|i58|) are descendant from the local zero modes on the topological charges 
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(instantons or dyons), and their density is therefore proportional to the den- 
sity of instantons (dyons). There are the instanton model [7,8] and the dyon 
model [9] of the QCD vacuum, which explain CSB in this way. 

Whether these models are realistic or not, is the open question, but the 
Banks-Casher relation (|i58|) holds independently of that, and if the method 



of the present paper proves CSB due to the field correlators (even in the Gaus- 
sian approximation), one should explain the origin of the quasizero modes in 



To do this we consider first the case of Abelian fields. As was stressed 
above in the paper, CSB is due to the correlator D{x), and the latter in 
the Abelian case can be connected to the correlator of magnetic monopole 
currents [1,2] 

< 3(3{x)35{y) > ^Qj.^QyJP^ partialxpdyp^^^^ ^^^^^ 

In the nonabelian case one can use the Abelian projection method (APM) 
[25], to separate in the field and field strength F^^ the monopole and 
photon part, and the part of "charged gluons" . The latter contributes around 
10% to the effective action. In this case one can connect the monopole current 
obtained by APM with the correlator D as in QlbUl ). 



Now for each magnetic monopole (or dyon) there is an infinite number 
of fermion zero modes, proportional to the length T of the world line of the 
monopole. Therefore the density of fermion zero modes per unit 3d volume 



and per unit of length along the world line is exactly as in ( |i58|) . On 



the other hand one can estimate this density from the 3d density of magnetic 
monopoles, which can be obtained from < j/3(a;, a:4)j^(0, 0:4) > integrating 
over Sx. (The correlator < j{x)j{0) > estimates probability of finding a 
monopole at the point x, if there is one at a: = 0. Integrating over d^x 
one finds the probability of having a monopole at a: = 0, while another is 
anywhere fixing X4 means that the probability refers to a given moment. We 
assume that one magnetic monopole yields one quasizero fermion mode per 
unit length of its world line - it is true for an isolated monopole, and we 
extrapolate this relation to the QCD vacuum as a whole.) Hence one gets an 
order of magnitude relation 
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J d^x < j(f , X4)j(0, X4) D{0)Tg (161) 
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where we have assumed for D{x) the form D{x) = D{0)f{^) and f{y) is an 
exponential or Gaussian with /(O) = 1. 

Finally, taking into account that D{0) < F(0)F(0) > one obtains 

< iIjiIj >« < F\0) > Tg (162) 

This estimate coincides with our result obtained from the quasiclassical cal- 
culation in the previous section. Numerically ( [TB2|) is - {SOOMeV)^, i.e. a 
reasonable order of magnitude. Thus the very existence of the "wrong" cor- 
relator D{x), violating Abelian Bianchi identity may bring about monopole 
currents and ascociated with those zero modes. 



8 The contribution of higher— order correlators 

The term in the cluster expansion of the effective action, proportional to the 
connected average of <C A^^{x^^'>)...A^^^{x^^'>) ^ contributes to the operator 
M the quantity 

iM^^\x^'\...x^^^)=^,,S{x^'\x^'^h,,...^,,^_,S^^^^^^^ 

^N^;LS^^'\-^^''^) (163) 

where we have defined 

«i='.,„K(i))...f',„^..K<"')» (164) 
and a(^4) = 1, Oi{^\^^) = %y, ^ = 1, 2, 3; A; = 1, ...n. 



One can identify in cumulant <C ... ^ the part similar to i.e. violating 
the Abelian Bianchi identity, namely for even n 

« ^.iMi(^^'^)-^.„.„(^^^^) »= n(^....w - 5,.,A.,.)D^''\^^'\ •••^^^^) 

i,k 

(165) 

Assuming for D the Gaussian form (the result which follows does not 
depend on that assumption modulo numerical factors) 

1 n 

D{u^'\ ...w(-)) = Doexpi--— E {u^'^ - u^'-'^f) (166) 



33 



one obtains at large x'''^\i = 1, ...n, \x^^^\ ^ Tg, l^*^*) — x^'^^ < Tg, z, A; = 1. 



.n. 



^™ = const|f I n(W-) (167) 



Now one can proceed in exactly the same way as we did for the bilocal 
Gaussian correlator, namely one performs Fourier transformation in 2:4"^ and 
for small P4^\p4'^Tg <C 1) one obtains the integral (we omit 7-matrices and 
Kronecker symbols for simplicity) 

,~(1) ,~(n-l) 
J ZTT ZTT 

n-1 

= iV(") n Ai(f^^\f^*+^^) (168) 

k=l 

where Ai(x, ^) is the same function as obtained in section 5, Eq.(81). Hence 
Ai is a moderated (5-function, which puts all x^^^ in the vicinity of each other, 
and M^"^) is the scalar function linearly growing with \x\ ~ Thus all 

conclusions about CSB and confinement of light quarks, which have been 
made based on the Gaussian correlator N^'^\ hold true also for any even 
order correlator N^'^^\ provided it contains the Z)-type term D^^\ 

Now the question arises: what is the contribution of all the sum 
Ej^i M{2k), and what is the result for the exact QCD vacuum? 

Here are two possible patterns. 

i) In the first case the higher order correlators are suppressed as compared 
to the Gaussian one. This may be true for the real QCD vacuum, and there 
are at least two evidences for the suppression of higher, n > 4, correlators. 
One fact concerns the behaviour of the static potential for static quarks in 
higher SU(N) representations. The gaussian correlator yields the quadratic 
Casimir operator as the coefficient on the potential, while higher correlators 
induce additional group structure as well [2] . The numerous lattice calcula- 
tion yield the firm evidence for the quadratic Casimir operator (see [2] for 
discussion and refs.), while other structures are not seen. 

Another fact is the study of the QCD string profile, i.e. the field distri- 
bution inside the string. This profile computed from the Gaussian correlator 
previously defined on the lattice [23] , coincides with the one, obtained by the 
lattice Monte-Carlo data [21]. Again no sign of the contribution of higher 
correlators. 
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Thus so far one can accept the hypothesis of the dominance of the Gaussian 
correlator, the corresponding picture is called the Gaussian stochastic vacuum 
[1,26]. 

ii) In the second pattern all higher correlators are important. One typical 
example, when this happens, is the dilute instanton gas model [7,8]. The 
effective Lagrangian looks like very similar to our starting expression, with 
the only exception: instead of averaging the product < A{l)...A{n) > over 
vacuum fields and expressing it through the field correlators as in one 
averages the instanton field A{i) over positions and color orientations of in- 
stantons. The result for the effective Lagrangian at large Nc was obtained in 
[27] 

n n _ 

(27r)4^(E^.)^r(n <fe))n^a.A-. + ^ ^ (169) 

i^i j=i 

Here is the standard instanton vector potential in the singular gauge, N 
is the number of instantons in the volume V4 and I ^ I implies a sum over 
antiinstantons. 

From (|1(39|) one easily finds the contribution of all terms with the product 
of n instanton fields (which is an equivalent to the n-th field correlator), this 
has been done in [27] with the result 

00 TV 

n n _ 

X (27r)4^(E9.)^r(n A\,{q,))^I^ I (170) 

In the sum ( |iYO|) all terms are important, and the sum can be computed 
explicitly [27] resulting in the equation, which was found earlier by another 
method [28]. 

M{p) = <p\{s- {A'r'r'ip > +/ ^ / (171) 

Here S is the averaged propagator 

S = \-Fr^ (172) 

p - iM{p) ^ ' 
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The nonlinear equation (|i7 i| ) can be studied in the low density limit, 

<C 1, A^c ^ oo, where p is the average size of instantons. One can argue, 
that the nonzero M{p), i.e. the CSB phenomenon, occurs for any density 
[27,28]. 

At the same time, as mentioned above, the confinement and hence the 
string between the light quark and the heavy antiquark here is missing, since 
instantons do not ensure confinement [2]. 

Thus the effect of higher correlators in the case of instanton gas model 
is twofold: on one hand higher correlators cancel in the sum for the string 
tension and hence destroy confinement, while it was present in the Gaus- 
sian approximation [1]. Remarkably this cancellation occurs only for integer 
topological charges as is the case for instantons, and is absent for magnetic 
monopole (dyonic) charges [29]. 

On another hand the higher correlators sum up in the new nonlinear equa- 
tion ( pn| ) which produces CSB as well as in the case of the purely Gaussian 
correlator, Eqs.(15),(16). 



9 Discussion and prospectives 

We have studied in this paper the simplest gauge-invariant system, contain- 
ing a light quark bound together with a heavy quark. This gauge-invariant 
setting allows for the appearance of the QCD string and the main subject of 
the paper is the influence of the string on the dynamics of one light quark 
at the end of this string. The vacuum correlator method has enabled us to 
introduce and to describe the string in the model-independent way, and the 
most part of the paper is denoted to the discussion of the bilocal Gaussian 
correlator. But this is because that simplest correlator already ensures the 
appearance of the string, and the most general correlators, discussed in sec- 
tion 8, do the same job, and as it was shown there, bring about the CSB in 
the same way as the Gaussian correlator. 

There have been used several simplifications. The first one is the large Nc 
limit, which allows to neglect additional quark loops (quenching approxima- 
tion) and to factorize averages of a product of loop integrals (e.g. in Appendix 
1). Correspondingly the problem effectively reduces to the one-flavour prob- 
lem, since the diagrams taken into account are the one-quark flow diagrams, 
(however with backward motion in time). This allows us to calculate the chi- 
ral condensate, which essentially constitutes a closed light quark loop, and by 
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choosing in S{x4 — ?/4, x, y) both times 0:4, 7/4 and coordinates x, ^ coincident, 
one disconnects the light quark loop from the heavy antiquark. 

So far the problem of several light quark flavours was not considered in 
the paper. As the next step of the same formalism one can introduce the qq 
Green's function and write for it the nonlinear equation similar to equations 
of the present paper, but taking into account both light quarks. Evidently 
the same formalism works for the baryon - 3q - state. This will be a subject 
of a future work. 

The main idea of the paper and of this future formalism is the selfinterac- 
tion of quark, which brings about nonlinear interaction kernel in the Dirac- 
type equation, and this nonlinearity breaks chiral symmetry. In this way the 
spontaneous symmetry breaking happens, and this seems to be a much more 
general phenomenon, than was assumed before. In particular, the same is 
the pattern of CSB in the NJL model. It is likely that the very concept of 
nonperturbative phenomena occurs in the same way - as the spontaneous 
breaking of scale invariance, and that happens via new solutions of nonlinear 
equations for gauge-invariant correlators. This is a fascinating way for the 
construction of the nonperturbative theory of gauge fields, in particular QCD 
[30]. 

Coming back to the topic of the present paper, the main task left un- 
completed is the exact calculation of the chiral condensate. The realistic 
estimates based on WKB have been given in section 6, but the logarithmi- 
cally divergent sum has been cut off in a reasonable, but approximate way. To 
calculate chiral condensate exactly, one should solve the nonlinear equations 
(15), (16) without local approximation, used in section 6 and valid actually 
at large \x\. This will be the subject of future publications. 
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Appendix 1 



Derivation of the effective Lagrangian (12) 

In this Appendix two ways of derivation of (12) are described: 

1) a simple way using the gauge (5), (6) 

2) a gauge-invariant derivation. In both cases we assume that is large. 
We start with the first one. Writing with fundamental color indices as 

Af = t^,A = l,...,N'^-l,a,b=l,...N, 

trt^t^ = ^Sab (Al.l) 

one has in the gauge (5), (6) 

Af: = l a{u)duF^^{u) (A1.2) 
Hence for the averaged product one obtains 

< Af{z)A^{w) >= titZ J a{u)du,j a{u')du', < F^^{u)F?,^,{v!) > (A1.3) 
Using now color neutrality of the vacuum one has 

/ < F^(u)F«(«') >= j^g' < F^{u)F'='{u') > (A1.4) 

The r.h.s. of the Eq. (A 1.4) can be expressed through D, defined in (10), 
namely 

< F^\u)F^{u') >= 2NcD{u - u'), (A1.5) 
since < trF{u)F{u') >= \ < F^{u)F^{u') >. Finally using the relation 

^ihl^cd = -^^adSbc - ^^abScd ^ -^^ad^c (Al.6) 

One obtains 

N r r 

< Af{z)Af{w) >= j^^SadSbcJ a{u)du J a{u')du'D{u - u') (A1.7) 
At large Nc Eq.(A 1.7) being inserted in Eq(4), reproduces (12). 
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2) In gauge-invariant version of relations (5), (6) one can insert there 
parallel transporters as follows 

^4) = / a^^\u)dui^^^' {z, z,; u, z,)Ff^'' {u, z^)^'^'\u, z^- z, z^) (A1.8) 

the same representation one can use for A^i{w). 

Now one can insert in both expressions (A 1.18) for A^j^{z^z^ and 
A^fiw^w^ the unity on both sides of F, namely 

24; 0, 24)^(0, 24; u, 24) = 1, for An{z,z/^) (A1.9) 

and 

1 = ^{w, W4; Ow;4)^(0, w^; Oz4)^(0, Z4; Ow;4)^(0, w/^; w, W4) (Al.lO) 

As a result one can write 

Af{z, Z4) = J a^^\u)dui^^''' {z, ZA] 0, ZA)Ff^^\u, z^- 0, z^)^^'\{), z^- z, z,) 

(Al.ll) 

A''^[w, w/^) = j a^'^\u')dUi,^^''' {w, W4; 0, 2;4)i^y'(w;, W4, 0, 2:4) $'^''^(0, Z4; w, w^) 

(A1.12) 

Now both F^ ^ and F^^ are referred to one point (0, Z4), namely we have 
defined 

F^'^'iu, Zi] OZ4) = ($(0, 24; 24)F(w, 24)$(m, 24; 0, 24))"'^' 

^^''^'(w;, W4; OZ4) = ($(0, Z4; w, W4)F{w, W4)^{w, W4; 0, z^))"''^' (A1.13) 

The average value of < A^A^' > (cf Eq. (7) of the main text) is expressed 
through the product of F in Eq. (A 1.13), and both F there are gauge 
tranformed as 

F(f , X4; 0, Z4) n^iO, Z4)F{x, x^; 0, ^4)^^(0, z^) (A1.14) 

and therefore the averaging of two F with this property produces the gauge 
invariant, namely 

< -2^4; Oz,)F^^^:iw, w^; 02:4) >= 

= ^ t < ^^(^(^' -^4; 02;4)F(w;, w^; Oz^)) >= (A1.15) 

c 

= — ^^^^ {Sii'^iitx' - Sif,'^i'n)D{u,w) 
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Therefore one obtains 

< Af{z,Zi)A';^{w,Wi) >= J a^^''\u')du'^, j a^^\u)duiX 

X Z4; Oz4)$"''^(0, Z4; w;4))$''^'(w;, w^- Oz4)$''^(0, 24; ^4) x (A1.16) 

1 

The insertion of this expression into Eq. (4) yields for the term bilinear in A 
^eff = jdzj dw{'4;^{z)j^'ijjbiz)){ijj+{w)-f^'^/jdiw))x 

{5,,^,'5u: - Si^>Si'^)Jlii' {z, w) (A1.17) 

Here we have defined 

= {^{Z, Z4; 0Z4)^{0, Z4, W, W4))ad (Al.18) 

Finally we note, that enters ^f, in the expression (A 1.17) which is 
assumed to be averaged over A^, and indeed, the averaged correlator (A 1.15) 
enters J-^i^ , which contains D{u — u'). Therefore the ^, in (A 1.17) are 
to be understood as averaged in the invariant combinations: 

< il;+{z)^ad{z, w)Mu^) >A, (A1.19) 

< Mz)'i'tbiz, ^)V^cH >A, (A1.20) 

Note also, that the invariant combinations (A 1.15), (A 1.19) and (A 1.20) 
can be averaged separately inside the common expression Cef / with accuracy 
of l/Ar2for large AT^. 

Finally, one can see that D{u, w) in (A 1.15) is defined in a gauge-invariant 
way, but depends not only on u, w, but also on the contour connecting those 
points. However one can see in the main text, that |2;4 — w;4| ^ Tg and in the 
limit Tg ^ one can replace the complicated H-shaped contour in (A 1.15) 
by a straight line connecting u and w, with an accuracy of Tg. For the 
straight-line contour the average in (A 1.15) and hence D{u,w), reduces to 
the function D{u — w). This is what is written in Eq.(lO) of the main text, 
and used throughout the paper. 
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Appendix 2 



The term A^^) and its role in the deconfinement transition. 



The term A^^^ defined in (10) iias tiie representation 

1 d 

^^^^ = o{^[(^ - u')i'5^^> -{u- + 



2 dui 

d_ 

du, 



+ ^^^[(^ - - - u')v5i^]}Di{u - u'), (A2.1) 

where derivatives act also on Di, so that the whole form (A 2.1) is a total 
derivative, and hence in the double integral (Pucfu' entering in the calculation 
of the Wilson loop [1] it contributes not to the area law term 
< W >^ exp{—a area), but rather to the perimeter term expCL where L is 
the length of the contour. 

The contribution of A*^^^ to the effective Lagrangian (12) is obtained re- 
placing {66 — 66)D in (13) by dS^\ given in (A 2.1). As a result one obtains 
J-^f consisting of a sum of 4 terms, which correspond to 4 successive terms 
in A(i), Eq. (A 2.1) 

=^ J^i^i' +^ J^^' +^ J^^i' +^ Jfifi'- (A2.2) 
We assume for Di the same Gaussian form as for D, namely 

Di{u) = L»i(0)exp(-^zV47fi) (A2.3) 

'Jm{z,w) = T^iDi{0){e ^ + 1) (A2.4) 

^Jkk' = a{u)du, a{u')du^,-^{{u - u')^'Di{u - u'))6kk' (A2.5) 



Jik = "^d^ Jo ^K^kt - Wkt')Di{zt - wt')- 

- i j^{zk - Wkt')Di{z - wt')t'dt'] (A2.6) 
Ju = - dui du'—[{u - u%D^iu - ^0] (A2.7) 
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'Jkki = 2 I c,iu)dui a{u')du'—[{u - u')k>Di{u - u')] (A2.8) 



^J^.^.'= J^'^.{z^w) (A2.9) 

We start now to discuss the contribution of different terms (A 2.4-A 2.9) 
to the kernel J{q,^) as in (20), and to the dimensionless kernels J^j^ as in 
(26-27). If one writes those in the form 

%M = ^^^^^(«7?)V;,'(Q,Q0, s = l,...,4 (A2.10) 

then the power u for each term ^ J would characterize its importance in the 
limit q; — ^ 0, e.g. terms with i/ = might contribute to the final equation 
(32), while the terms with i/ > are suppressed and vanish in the limit 

3 
2 



We start with ^ J44. The first three terms on the r.h.s. of (A 2.4) yield v = 



and therefore are not important in the limit a — ^ (i.e. at large distances). 
The last term on the r.h.s. of (A 2.4) is constant, and a simple calculation 
yields for that u = 0, therefore one should look at the effect of this term 
more closely. We shall show now that the constant term in J cannot yield 
the selfconsistent nonzero solution of equations (20), (21) or (31), (32). To this 
end we insert the constant J{z,w) or equivalently J{q,^) = const6{q}6{^) 
on the r.h.s. of (20), and realize that solutions for M, S should have the form 

S{PA,P,0) = a{pA,p)6{P + ^){27Tf (A2.11) 
From (20), (21) one then gets an equation 

^ip4,p} = const^=^=== (A2.12) 

/p^ + II\P4,P) 



which has only the trivial solution fi = 0. 

For (A 2.5) one gets vi^Jik) = 1, this term is unimportant. The term ^744, 
Eq.(A 2.7) contains the full time derivative, which being Fourier transformed 
as in (32) yields additional factor (P4 — P'^'^aTg which can be neglected. For 
the rest kernels, Eqs. (A 2.8) and (A 2.9), the corresponding v = 1. 

Thus the analysis of the nonconfining correlator Di, brings one to the 
conclusion that it cannot yield a selfconsistent solution for the scalar kernel 
M, and hence cannot break chiral symmetry. 
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Appendix 3 



Properties of the kernel Jik 



As defined in (13), Jik is 

Jik = Jq duia^iu) a^iu)dUkD{u - u) = 

= ZiWk dtanit) dt' a^i{t')D[zt — wt\ Z4 — W4) (A3.1) 
For tlie Gaussian form of D, one has 

D{zt - wt', z, - w,) = DiO)expi- )exv{^- ~ f ^^' ) (A3.2) 

g g 

The last factor can be rewritten as 

(A3.3) 

The integral / dtdt'r]{z, w) has different asymptotics in 3 different regions 
of the z, w space. 

3)\z-w\ ~ \z + w\ ~>Tg 
In the region 1) one has 

4- = z,w,D((^)e ^1^1^(1 + 0(-^)) (A3.4) 

Jik — '^Jik (A3. 5) 
In the region 2), 3) it is convenient to rewrite the integral Jf^ as 

4 = z,w,D{0)e-'^ j^j^ di f \ dXe-^'-i'-'^'''^ (A3.6) 

12'. A- till It: — in\ JO J-ct 



\z + w\, 


1 — * 

\z — 


w\ 


\z + w\, 


1 — * 

\z — 


w\ 



where we have defined 



(Z^ — W'^) \Z + W\ ^ \Z — W\ / . o r,N 
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For c ^ 1 and ^ 1 one can expand the integration region in A to 
(—00,00) and shift the variable as A + ^cosx = A', which yields for the last 
two integrals in (A 3.6) 

[ V^^o, ^o\sinx\ ^ 1 

Insertion of (A 3.8) into (A 3.6) finally yields in region (2) 

4^ = z,WkDiO)e i (A3.9) 

For i^o|s^^xl ~ 1 one using (A 2.8) recovers the result (A 3.4). 

In the region 3) the integrals in (A 3.6) are of the order 0(1) unless cosx = 
±1, and the estimate of Jf^ is typically therefore 

Jf^ = z,WkD{Q)e "^ .^^ ^|^^_ ^J {\{z + w\/\z- w\) (A3.10) 
where we have defined 

/ = ^ /° Jll^ dAe-^'-^'-2^^-^^ (A3.11) 

For cosx = when z and w are (anti) parallel and both large in 

modulus, one can choose coordinate (r, r') in (A 3.2) such that r = t — ^t', r' 
orthogonal to r, and integration over dr yields the factor as in the case 

Finally one can write for J-^, noting that 

W = 4T2(-i^ - -^) (A3.12) 

therefore since in the integral (A 3.6) both A and ^ are of the order of one, 
hence one obtains in the regions 2), 3), that 

'^ik ~ (A3.13) 
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Appendix 4 



Corrections to the kernel (48) at finite Tg 



In the calculation of M^p^, z,w), Eq. (44), one should do the integral 

Nip,) = / ^ (A 4.1) 

In section 4 the limit Tg ^ was used in which case the answer was 
given in (47). The limit however depends on the boundary conditions at 
infinity, and therefore should be done with care. Moreover we calculate in 
this appendix corrections due to finite values of Tg. To proceed one can 
represent N{p/^ as follows 



N(P.) = r ^^"'" Kft^) (A 4.2) 

where 

^^^'^ = J 27r(y4 - isn) ~ ^ -^^(-/.4)^(-^n)e-l^"^^l, < ^-^^ 
Inserting (A 4.3) into (A 4.2) one gets 

Nip,) = f dh,e-i~''^''''"''' (A 4.4) 

Taking into account, that 

r dhe-^-"' ( '"''^P'''^,] = (A 4.5) 

^0 \^ sm{p4h) J ^ ' 

1 \ / ^^Vikl - ^P4)^T2)[1 - $((7 - zp4)T,)] 

-i j 4 1 ±exp((7 + ip4)%')[l - $((7 + ^P4)T,)] 
one obtains the following expansion of N at small T^, 

JV(p4) = ^i^f<^=>(l-^ + (4-rf)7;^-^(.«Sn£„)r,+0(T,^) (A 4.6) 



In the opposite limit, when \pA^g ^ 1 or/and Iffnl^^ ^ 1 one has from 
(A 4.5) 
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Appendix 5 



The limit of the large mass m 



In the equation for S 



{-id-im)S{h4, z, w)-J iM{h'^, z, z)S{h4-h'^, z, w)dh'^d^^ = S{h4)S^^\z-w) 

(A5.1) 

for large m ^ M one can neglect in the first approximation M and find the 
free Green's function 5*0 for the massive quark, 

S'o(/i4, z, w) = {im — id)f{z, w), (A5.2) 

Since we shall treat differently = Z4 — W4 and z — w, another representation 
is useful, 

5o(/i4, z,w) = - --^^_^{e-^^^/^^(/l4)(m + 74^^^^ - iff) 

+ e^'^/^+^'e{-h^){m - 74^^^ + f - ip^)} (A5.4) 
The form (A 5.4) allows the expansion in powers of 1/m, which yields: 

So(h4,z,w) = SQ^\h4, z,w) -\ 5o^V/i4, £, w;) + ... (A5.5) 

m 

where 



Si'\h^,z,w) = ^^('H^-^ZJ){e-'^^"^^(/i4)(l+74)+e'^^"^^(-/i4)(l-74)} (A5.6) 
Insertion of (A 5.6) into the expression for M 

iM{h4, z, w) = J£(/i4, z, w){j^S{h4Z, w)jf,Sik - ikSii) (A5.7) 
i.e. replacing S by yields 

zM(/i4, £, w) = '-J[z, £)^(3)(£_^){^(/i4)e— ^^[2-i(l-74)]+^(-/i4)e-'^i2-l(l+o 

(A5.8) 
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where we have defined 

J{z, z) = z^ dtdt'D{z{t - t'), /14) (A5.9) 
For the Gaussian form of D one has at large \z\: 

J{z,z} = ^e~^ (A5.10) 
Insertion of (A 5.8) into (A 5.1) yields 

(j\z\ f 

x{e{h',)e-^^'^[2 - i(l - 74)] + e{-h!,)e^'''^[2 - ^(1 + -i,)]]S{h, - h',z, w) = 

= S(h4)6^^\z-w) (A5.11) 

To get rid of the /i4 dependence, we multiply both sides of (A 5.11) by e^^^ 
and integrate over dh^. Defining 

S(z,w) = r e'^^'S(h4,z,w)dh^, (A5.12) 
J— 00 

one obtains from (A 5.11) an equation for S, 

{-ijd - im{l - j4))S{z,w) - ia\z\{[l - hi -74)]+ 

D 

+ [1 - i(l + 74)]J(2mr^)}5(£,w;) = 6^^\z- w) (A5.13) 
where we have defined 

^ X > 1 



Jix) = expix')[l - $(x)] = ^ I [ (A5.14) 

where ^{x) is the error function. 

In arriving to (A 5.13) the change of variables h4^ — h'^ = u, and integration 
over dh^du helps to factorize out S under the integral. 

Till now we only have made an expansion in 1/m keeping the leading 0(m) 
and next-to-leading 0{mP) terms, in particular the kernel M ~ (7\z\ is 0(m°) 
when one inserts in (A 5.7) So instead of S, implying that (j\z\ <C m. 
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Therefore the product rriTg is arbitrary and one can consider two hmiting 
cases as in (A 5.14). For mTg ^ 1 one can neglect the term with J in (A 5.13) 

and obtains an equation for S = 

-adS-+ - ia\z\S++ = 5^^\z - w) 

-ia\z\S+- - adS— = 

- 2i ^ - 
adS++ — 2imS-+ — — (7|i'|5_+ = 

o 

adS+- - 2imS— - i^cr|£|5__ = 6^^\z- w) (A5.15) 
To the leading order one obtains as expected 

-2 

(-^ + (t\z\)S++ = i5^^\z- w) (A5.16) 
2m 

In the opposite limit, when mTg <C 1, one has from (A 5.13) 

r 

[-ijd-imil - 74)]5 - i-a\z\S = S^^\z- w) (A5.17) 

o 

and in leading order of 1/m one obtains instead of (A 5.16) 

+ l^\^)S++ = iS^'\z- w) (A5.18) 

We conclude this appendix by discussion of the role of magnetic and elec- 
tric correlators ^ on the properties of confinement of the light quark. 

Having in mind that in (A 5.7) the term 74574 enters with ~ D^, while 
7j57i is multiplied with ~ D^, one can see that for = 0, 7^ one 
would have in (A 5.8) 

iM^ = '-J^iz, z)S^'\z- w){eih,)e-^^^^il - 74) + ei-h,)e^^^'^il + 74)} 

(A5.19) 

To the leading order in 1/m one obtains the free equation instead of (A 5.16) 
implying that magnetic correlator does not lead to confinement. 
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